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Abstract 

We establish the analytic non-integrability of the nonholonomic ellipsoidal rattleback 
model for a large class of parameter values. Our approach is based on the study of 
the monodromy group of the normal variational equations around a particular orbit. 
The embedding of the equations of the heavy rigid body into the rattleback model is 
discussed. 
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1 Introduction 

The rattleback's amazing mechanical behaviour can be described as follows: when spun on a 
flat horizontal surface in the clockwise direction this top continues to spin in the same direc- 
tion until it consumes the initial spin energy; when, however, spun in the counterclockwise 
direction, the spinning soon ceases, the body briefly oscillates, and then reverses its spin 
direction and thus spins in the clockwise direction until the energy is consumed. The first 
mathematical model of this phenomena belongs to Walker (1896) who studied the linearized 
equations of motion and concluded that the completely stable motion is possible in only 
one (say clockwise) spin direction. This classical explanation of the rattleback's behavior is 
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incomplete since it does not reflect the global dynamical effects explaining the transfer of 
trajectories from the vicinity of the unstable solution to the stable one. To analyze thor- 
oughly this question we propose to study the nonholonomic equations of the rattleback in 
the complex domain and particularly to study the existence of additional analytic first inte- 
grals. It has been observed in many mechanical systems that non existence of analytic first 
integrals is usually associated to complicated chaotic behavior of trajectories of the system. 
We mention that actually only numerical evidence for chaos in the rattleback systems has 
been observed [I]. 

In our case the rattleback represents a full ellipsoidal body whose center of mass P coincides 
with its geometric one. All vectors are defined in the body fixed axes (with origin in P) 
coinciding the with principal geometric axes of the ellipsoid 
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BM = % + % + % = l. (1.1) 

Let ( , ) denote the euclidean scalar product in M 3 and let 1 1 • 1 1 be the corresponding norm. 
In addition, let [, ] denote the cross product in M 3 . 

As shown in |3] , the nonholonomic equations of motion can be written in the following form 



m [r, [lu, [lu, r])) + mg [r, 7] - m [r, [lu, —)) , (1.2a) 

(1.2b) 
where 
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is the inertia matrix; m - mass of the body; g - gravitational constant; lu = (lui, u> 2 , uJs) T - 
angular velocity; 7 = (7i,72,7s) T - unit vector normal to the body's surface at the point of 
contact; r = (ri,r 2 ,r 3 ) T - position of contact point. 

To observe the rattleback behavior, one has to impose the following condition on G (see 
e.g. pQ) and bi 

£11 S 12 
£12 S22 

E 33 = I 3 , h^b 2 , 

where I\, I 2 , I3 are the principal components of the inertia tensor whose principal horizontal 
axes are rotated by the angle S > with respect to the ri-r 2 -axes of the ellipsoid. 



_ duj , ,du 

d'j 
~dt 



—[lo, Qlo) 
[7, w] , 



Ii cos 2 5 + I2 sin 2 5 (ii — I 2 ) cos 5 sin 5 
(1 1 — I2) cos S sin 5 I\ sin 2 5 + I2 cos 2 S 



f lA) 
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Unless otherwise mentioned below we will assume that 

£y>0, £3,3 < E M + s 2, 2 , 61 > & 2 > 0, 63 >0, m>0, g>0. (1.5) 
Using (jl.lj) . the vector r = #(7) is defined by solving 7 = — VF/||VF|| for r, which gives 

r« = Hi) = -f= 2 2 ~ b ; 2 T 2 ,22 ' * = 1 > 2 ,3. (1-6) 

Solving (jl.2aj) for Co and using $Oj) we transform the nonholonomic rattleback equations to 
the standard form 

(1.7) 

c?7 

where F = (Fi, F 2 , F 3 ) T is a vector field rational in the variables u, 7 and s = v^iTi + b\^\ + 6373. 
Let (a;, 7, s) G C 7 . Then the vector field (|1.7|) . as a function of (cj, 7, s), is analytic in the 
domain D = C 7 \ (5 U £) where £ = {(a;, 7, s) G C 7 : s = 0} and S C C 7 is the surface on 
which the determinant of the matrix U given below is zero 

U = — m[r, [r, ■]] = © + m((r, r))Id — r <g> rj . (1-8) 

This matrix appears when we solve (jl.2a|) to find Co. For an arbitrary non-zero vector m6 1 3 

(u, Uu) = (u, Qu) + m(u, u) (r, r) — m(u, r) 2 
= (u, Qu) +m\\u x r|| 2 > 0, 

so that U is positive definite and hence det(U) = never occurs in the mechanical case. 
The system 11.7)1 always possesses two first integrals (see |3]): 

777 I 

H = -||[^,F(7)]|| 2 + -(u,Qu) - mgiR^lj) = h, h G R, (1.9) 

- the energy; 

G=(7,7)=7i 2 + 7 2 2 + 7 3 2 = ^ (1.10) 

- the geometric integral. 

We see that H(u,j, s) (after introduction of s with the help of p. 6)1 ) and £(7) are analytic 
functions of u, 7, s in T>. The natural question arises whether the equations (jl.7)l can have a 
third first integral analytic (meromorphic) in (a;, 7, s) G C 7 and functionally independent of 
H and G. We remark that the absence of first integrals in the meromorphic case is generally 
harder to prove than in the analytic one. 

The paper is organized as follows. In Chapter El we describe one particular solution of the 
Rattleback problem. Chapter El discusses some properties of the normal variational equations 
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and the monodromy group generators of this solution. Chapter 0] examines the case in which 
the normal variations equations have logarithmic branching points. Under these conditions 
we show the non-existence of additional meromorphic first integrals (Theorem |2J). Finally, 
Chapter El contains the proof of our main theorem about the analytic non-integrability: 

Theorem 1. Assume that the conditions ()1.5|) and ()5.3|) hold. Then the energy H and the 
geometric integral G are the only integrals of the rattleback problem (jl.7j) which are analytic 
in (u>, 7, s). 

Rigid body limiting case 

The rattleback equations ()1.2aj) - ([1.2b[) formally contain the equations of the heavy rigid body 
in the singular limiting case m — > 0, mg -/-> 0. Moreover the functions Riiri) are constants 
denoted by r$ again, now designating the position of the center of mass in the body frame. 
By a rotation the tensor of inertia G can be diagonalised so that S 12 can be set to zero 
in this case, thus violating (jl.5j) . Under these assumptions the system is Hamiltonian with 
another well known integral 

L = (7, Qu) = 1x<JJ\Ii + 72^2^2 + 73^3 • (1.11) 



2 The invariant manifold 

We consider the vector field (jl.7j) as a function of six variables (00, j). Let I > (in the 
mechanical case we can set 1 = 1). For (u, 7) G M 6 the square root in s is assumed to 
be always positive. It is straightforward to check that the system (jl.7j) has the invariant 
manifold M = {(a;, 7) G M 6 : uj\ = u> 2 = 0, 73 = 0}. The mechanical sense of the motion 
on M is quite clear: it corresponds to rolling (or oscillating) of the body on the line of 
intersection of its surface with r 3 = plane. We note that this invariant manifold exists 
because of the assumption that the r 3 ellipsoidal axis coincides with one of the principal 
inertia body axes. 

Complexifying (jl.7j) we denote by Ui t h C M C C 6 its orbit in M corresponding to fixed 
(I, h) G C 2 defined as the intersection of the two algebraic surfaces 

f m ^2\ + i¥ } + + 2 mg JbU + bhl = 2heR, (2.1a) 

V l7l + °272 / 

7l 2 + 72 2 = /GM. (2.1b) 

We assume now that h 7^ and I = 0. Obviously the solution obtained in this way cannot 
be real; it is non-mechanical. 

Let T C C 6 = (^,7) be the complex curve defined by 

p 2 — a 2 p 2 — a 2 

^3=P, 71 = — g — ,l2 = i — Tj — ,Wi = w 2 = 73 = 0,l)GC\{-a,+a} 1 (2.2) 



6o> 3 2 
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where 

a 2 = m 2 " v , p = - ^ffi^ -, ^ C > (2-3) 

This is the parametrization of Uo,h- 
Let Po = (cuo, 70) G T such that 70 7^ 0. 

The function s = a/^i7i + ^272 + once the square-root branching is fixed, is analytic in 
a small neighborhood Up Q C C 6 of Po- We can analytically continue s and all its derivatives 

p 2 — a 2 

along T with the help of the parametrization above: s((u,j) G T) = s(p) = a — where 

= \Ai > 0- m particular, it shows that s is single-valued on T and hence is analytic 
in a small neighborhood By C C 6 of it. Thus, the vector field (jl.7j) restricted to By is an 
analytic function of the complex variables (^,7) and has the invariant curve V. 

We use the reparametrisation bf = Pi(pi — P2), with pi > p 2 , Pi = const > 0, which makes 
(T — px — P2 polynomial on the solutions we are considering. In particular for this solution 
r = (— pi, — ip2, 0) T , which is a constant. 

Since in the rigid body limiting case r is not proportional to 7 the invariant manifold M only 
exists under the additional assumption that r 3 = 0. 



3 The variational equations and its monodromy group 



The relation between the parameters p and t can be easily found by substitution of the 
parametrization ()2.2|) in one of the equations (jl.7|) . This gives 

v 2 — a 2 

dp = i P dt. (3.1) 

Remark 1 . The change of time (j3.1|) has infinitely many sheets. Nevertheless, one can always 
replace the vector field F(X), X = (uj^) t defined in (JUj) by F = (F/(2~H(p 2 - a 2 )), 1) T 
i.e. consider the new autonomous system of differential equations 

dX F(X) dp = 

dr 2-H(p 2 -a 2 y dr ' 1 ' 

which will have the same particular solution V defined by ()2.2|) where p is replaced by r. 

Obviously, an autonomous analytic first integral of (jl.7)l gives a first integral of the same 
type for the vector field (|3.2jl . One sees also that F is analytic in the neighborhood of V. 
The further analysis based on the variational equations of (jl.7j) or ()3.2j) will be essentially 
the same. 

Let Ui, 7j, i = 1,2,3 be the variations of variables Ui, ji, respectively. We use the following 
notation for coordinates of the variation vector V = (X, Y) T = ((£>i, £>2> 73), (71, 72, ^3)) T - 
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In these variables, the variational equations of (J1.7)) along T, after the substitution (J3.1)) . 
have the block diagonal form 



dX 
dp 



M O 
O N 



X, X G T r B r . (3.3) 



where M(p), N(p) e GL(3, C(p)) and O is the zero 3x3 matrix. 

One can show that the derivatives of the first integrals H and G with respect to the variables 
^1, ^2) 73 vanishes along T. So, the linear first integrals H = (dH(T), Y) and G = (dG, Y) 
are not useful in solving the first block system 

| = MI, (3.4) 

dY 

whereas the second subsystem — — = NY can be completely solved in radicals with the help 

dp 

of H and G. 

As seen from the parametrization (J2.2)) . the equations ()3.4|) are the normal variational equa- 
tions (see [Zj for definition) of the system (jl.7j) along the orbit V . The following lemma 
holds. 

Lemma 1 (|6J,|7J). Let us suppose that the system (jl.7|) has a third first integral H^u^) 
which is analytic (meromorphic) in the neighborhood By and functionally independent of H 
and G. Then the monodromy and the differential Galois groups of the normal variational 
equations ()3.4|) have a non trivial polynomial (rational) invariant. 

We can put the linear system (J3.4)) into the Fuchsian form with help of the rational trans- 
formation 

X = diag(p, p,p 2 — a 2 )x , (3.5) 
where x = (xi,X2,X3) T G C 3 . The new system takes the following form 

dx ( A A B\ 

— = T(p)x=( + + — )x, peTZ = C\{-a,+a,0,oo}. (3.6) 

dp \p — ap + apj 

Clearly, the monodromy groups of the systems ()3.4|) and ()3.fi|) are equivalent. The constant 
3x3 matrices A, B are given in Appendix A. 

Remark 2. One sees that the above linear equations are invariant under the transformation 
p i — > —p. That in particular allows to reduce the number of finite singularities up to two via 
introducing the new time p = r 2 . Nevertheless, we prefer to keep this symmetry inside the 
equations and to use it later in the study of the monodromy group of ()3.6|) . 

Proposition 1. The general solution of ()3.6|) is meromorphic in a neighborhood of p = 0. 
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The proof follows from the fact that the matrix M(p) in (|3.4j) is holomorphic in the neigh- 
borhood of p = as seen from (|2.2|) . ()1.2aj) . In particualr, we have 

Spectr(S) = {0,-1,-1}. (3.7) 

That can be verified directly with help of formulas for B given in Appendix A. 

We will need some technical results. The following proposition will play a crucial role for 

our non-integrability results below. 

Proposition 2. The characteristic polynomials P(x), Poo(x) of the residue matrices A at 
p = ±a and = —2A — B at p = oo always have non-real roots. 

Proof. A direct calculation shows that P(x) is of the form 

p(x) = x 3 + x 2 + e x x + O , (3.8) 

where with the definitions 

*■ = £ ' * d = *■ " (E " + E * ~ E33,Icl " = ft) = (-£) ' (3 ' 9) 

the coefficients of P are given by 

9 d = det^ d + m(v,^ d v) (3.10a) 

6 d B x = det # n + m{(v, V n v) + p 3 (£nPi - £22/02 + iSi 2 (pi + to))} (3.10b) 

#A = 0A - mp 3 (pi - pa) (£n + S 22 - £ 33 ) . (3.10c) 

Therefore 

Im-^— = —, 2Sl2PlP2 -, (3.11) 

0o-0i p 3 (pi - p 2 )(£n + £ 22 - X 33 ) 

which is obviously non-zero according to (|1.5jl . A direct calculation shows that Poo{x) is of 
the form 

Poo (x) = x 3 - 4x 2 + Xix + xo (3.12) 
In the notation of the previous proposition and the Appendix the linear coefficient is given 

by 

Xl = 5 + Aa\ + 4a 2 a 3 + 4(c 3 - (3.13) 
The coefficients of P(x) and Poo(x) are related by 

12(0 O - 0i) + *o + 3xi = -9. (3.14) 

Therefore 

Im = (3.15) 

Xo + 3xi + 9 6p 3 (pi - p 2 )(Sn + S 22 - E 33 ) 

□ 
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In the rigid body case the matrices A and B become 



/ i(E 33 - S 22 )/E n i(E u - E 33 )/E 22 

(—1 — 2imgT 2 /En 
-1 2im#n/E 22 


Therefore 

Spectr(^4) - ' / (^n — ^33X^22 — £33) / (^11 — S 33 )(S 22 — S 33 ) 




^11^22 V ^11^22 



I 



For the Kovalesvkaya case En = E 22 , E 33 = 2En the eigenvalues become {—1, 1, 1} and 
{—2, 3, 3} at infinity with a non-trivial Jordan block. 

The following proposition shows that the system ()3.6|) is not solvable in the Lappo-Danilevsky 
sense (see j2]). 

Proposition 3. The residue matrices A and B do not have common eigenvectors. In par- 
ticular, they do not commute. 

Proof. The proof is a straightforward, but since the proposition will play an important role 
later, we indicate some basic steps of it. A direct computation gives for the eigenvectors Vj 
and the corresponding eigenvalues A» of B: 

Vi = (l,vi,v 2 ) T , Ax = 0, 
V 2 = (l,0,0f, Ai = -1, 
V 3 = (0,l,0f, A 1 = -l, 

where Vi are some known expressions. 

One verifies that the matrix R = (fi,f 2 ,t> 3 ) is not singular once the condition (jl.5j) holds. 
Under the same conditions it is easy to show that V 2 and V 3 are not eigenvectors of A. Here 
is one method to prove the same statement for V\. We consider the conjugation A = R~ 1 AR 
of A by R. It is sufficient to show then that 7^ or Ai )3 7^ 0. That can be done quite 
easily once the condition (|1.5J) is fixed. Otherwise, if C = [A, B], then the non-commutativity 
property follows directly from C 3 i — 1/(3^0. 

□ 



Fixing a basepoint e £ 1Z, one defines the monodromy group G of the system ()3.6|) as the 
image of the fundamental group tt(TZ, e) given by the analytical continuations of the funda- 
mental matrix solution E(p), E(e) = Id of (|3.fij) along all closed paths T G ir(lZ,e) starting 
from e. One verifies that tr(A),tr(B) e Z. Hence, G C SL(3, C). According to Proposition 
the group G = (M + , M_) is generated by the local monodromy transformations M + , M_ 
around the singularities p = a and p = —a, respectively. 
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4 On the non-existence of additional meromorphic first 
integrals in the logarithmic branching case 



The next proposition will show that the logarithmic branching of solutions of the normal 
variational equations (J3.6j) is not compatible with existence of a new meromorphic first 
integral of the rattleback problem (jl.7|) . 

Theorem 2. Let us assume that the general solution of the variational equations (jH.fi j) has 
logarithmic branching around one of the singularities p = ±ct. Then the rattleback equations 
do not have any new meromorphic first integrals. 

Proof. By a suitable linear transormation of the general solution of (j3.6j) . we can always 
reduce one of the monodromy matrices around p = a or p = —a to its Jordan form. In 
the case that the general solution has logarithmic branching points, at least one of the 
monodromy matrices M+, M_ has a non-trivial Jordan block. 

Since p = ±a enter in (jH.fi)) in a symmetric way, it is sufficient to consider the case of p = a. 
Firstly, we assume that M + is of the form 



Mi 



q 1 
q 1 
q 



(4-1) 



Since M+ e SL(3, C), one has q 3 = 1. 
Let 

Spectr(A) = (Ai, A 2 , A 3 ) 
As follows from (j3.6j) (see e.g. [I]) 



Spectr(M+) = Spectr(M_) = {e 2 ™ Xl , e 2 " A2 , e 2wiX3 } . 



(4.2) 
(4.3) 



Thus, in the case (j4.1j) all eigenvalues of A must be real (even rational) numbers. This is 
impossible according to Proposition 

We suppose now that M + is of the form 





u 


1 





M+ = 





u 













u~ 2 



(4.4) 



In this case we prove the following lemma. 

Lemma 2. Let R(x), x = (xi,X2,xz) T be a rational invariant of G. Then R = R(x2,xs). 
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Proof. We write R(x), degP = m as follows 

i 

R = ^r , I, k G {0,1,2,...}, m = l-k, (4.5) 

E »3~ J ^j( ar i> a; 2) 
i=o 

where Pi(xi,x 2 ), Qj(xi,x 2 ) are homogeneous polynomials of degrees i and j, respectively. 
Let us assume that at least one of the polynomials Pj or Qj depends on x\. So, for example, 
we will find < p, s < I such that the term x^P s (xi, x 2 ), dP/dx\ ^ 0, deg(P s ) = s is a 
semi-invariant of M + , i.e. it satisfies 



MlP s (xi,x 2 ) = c n P s (x 1 ,x 2 ), ceC*, WneN, M 4 



u 1 
it 



(4.6) 



Since P s (xi,x 2 ) is homogeneous of degree s, we can put it in the form 

P s (xi,x 2 ) = a x{ + aix\~ 1 x 2 + • • • + a s x s 2 . (4.7) 

We observe that M™(xi,x 2 ) T = (u n Xi + nu n ~ l x 2 , u n x 2 ) T . Consequently, the polynomial 
M™P s (xx,x 2 ) will have its coefficient of x| equal to 9 n = J2l=o ct r n s ~ T u 3 ^ 1-1 ^' . Let q be 
the smallest index such that a q ^ 0. Then the coefficient before x\x s 2 ~ q in M™ P s (x\,x 2 ) is 
rj n = u nq a q . Clearly, as n — > oo, the asymptotic behaviors of 9 n and 7? n are different. So, the 
equality (|4.6|) cannot be true. That finishes the proof of Lemma [21 □ 

Lemma 3. All rational homogeneous invariants of M + are functions of the invariant x\x^. 

The proof follows easily from the condition that Spectr(A) is non-real. 

Lemma 4. Let us assume that the monodromy group G = (M_,M + ) has a rational homo- 
geneous invariant R(x). Then x\x% is the invariant of G and R = {x\x^) m for a certain 
m G Z*. 

Proof. Let R(xi, x 2 ,x 3 ) be a rational homogeneous invariant of G written as 

u 'Jl,l2,«3 -^l ^2 A 3 

= h+h+is=M {0,1,2,...}, (4.8) 

°\n J2 J3 x l x 2 x 3 



where M,N are non- negative integer numbers and all coefficients Oj, 6, are different from 
zero. 

Firstly, we assume that R reduces to the single term 



R = a%x$x l i, keZ, (k,l 2 ,l 3 )^ (0,0,0). (4.9) 
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Then li — according to Lemma 121 and, as follows from Lemma E3 there exists m G Z* such 
that 

R=(x 2 2 x 3 ) m , (4.10) 

that in turn proves the result. 

If R is not of the form ()4.9|) . it can be written as follows 



R 



ixjy ^ t/y 2 **'3 1 2 3 ~~ 

„^Pl ™P2 J3 , ... 

to/ j x 2 x 3 t - 



ri,qj,Pk e {0, 1,2, 



} 



(4.11) 



where a, 6, c 7^ and (7*1, r 2 , r 3 ) 7^ 92, 93), S r i = Qi = M - 
We note that, as seen from (|4.11jh the terms X -y X<2 *^ 

q 3 and 

^ 1 ^2 ^3 £11*6 multiplied by the 
same constant under the action of M + . Therefore, the division by shows that 

r(x) = Xi 1 x 2 2 x 3 3 , ki = Ti — q i} (ki,k 2 ,k 3 ) 7^ (0,0,0), ^ fcj = is the invariant of M + . We 
have ki — according to LemmaEland hence r(x) = (£ 2 /£3) fc2 , fc 2 7^ 0. This is impossible 
according to Lemma El 

Thus, we can assume R to have the form (|4.10j) . This immediately fixes the monodromy 
transformation M_, also preserving R, as follows 



M_ 



U\ n k 
u 2 
u 3 



U1U2U3 = 1, Ui,n,k G C . 



(4.12) 



We will show that iti = w 2 = m and 1*3 
Indeed, since Spectr(M + ) = Spectr(M_ 



u 



{u, u, u 2 }, we have necessarily u\ = u or u 2 



it. 



Let Ui = u, u 2 = u 
unity. Let now u\ = u 



u 3 = u. With help of R we obtain w 3m = 1 and so u is the root of 
U2 = u, us = u. Then u 3m = 1 and we conclude as above. Hence, 



the monodromy transformation M_ takes the form 





u 


n 


k 




M_ = 





u 





(4.13) 










u- 2 





The rational function R = x\x 3 is clearly invariant under the action of M_ and M + . 



□ 



We shall follow the approach close to that of Tannakian (see e.g. [S]). 

Let E(p) be the fundamental matrix solution of the equations ()3.6|) and S _1 (p) = (Si, S 2 
where £, are linearly independent vector functions. It is known (see e.g. p], p. 246) that if 
R(x) — x^x^ is db polynomial invariant of the monodromy group G then 



S 3 ) T 



i{ P ,x) =R(E( P )- 1 x) = (j: 2 ( p ),x) 2 (j: 3 ( p ),x) 



(4.14) 
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will be a first integral of (|3.6j) invariant under the action of G, i.e. single-valued as a function 
of p G 1Z. One can express this fact by stating that all coefficients aij^ip) in the expression 
for I, 

are rational functions of p. 

Since the polynomial / given by (j4.15|) can be factorized as ()4.14|) . there exist two 3- vector 
functions A(p), £>(p), algebraically dependent on p, such that / becomes 

I(p,x) = (A(p),x) 2 (B(p),x). (4.16) 

The structure of the monodromy group generated by the transformations (J4.4j) . (J4.13|) sug- 
gests that the system (|3.6|) has two linearly independent particular solutions X 12 (p) of the 
following form 

Xi(p) = (p- a) Xl (p + a) Xl Riip), e 2mXl = u, Ai G Spectr(A) , (4.17) 

X M = (p - a» + «)*«,(,,) + u-%W + n!^|±^) , (4.18) 

where Ri^(p) are vector functions rationally dependent on p. 
Lemma 5. VFe /wwe (A, Ri) = = 0. 

Proof. Plugging the solution (|4.17|) into the integral ()4.16|) we obtain 

I(p,X x ) = (p-ai) 3Al (p + a) 3Al (^ ;J Ri) 2 (i3, R x ) = c = const . 

In the case c ^ 0, since A, B, R± are algebraic with respect to p, the last equality implies 
Ai G Q. According to (|Q)l 

Spectr(A) = {A, A + fc, -2A - fc - 1} , (4.19) 
for certain A G C and k G Z. 

So, it is clear that if Ai G R, then Spectr(A) has to be real. According to Proposition |2] we 
therefore get (A, Ri) 2 (B, R x ) = for all p. 

We shall consider the case (A,Ri) = and (B,Ri) ^ 0. Obviously, the shift I(p,x) = 
I(p, x + Xiip)) of the first integral I is again a first integral of the system (|3.6|) and 

J = J + (p - a ) A i (p + «) A i (£, (A x) 2 = const . (4.20) 

Since the equations ()3.6|) are homogeneous, with help of ()4.2(Jj) . ()4.16|) . we derive the first 
integral 

(p - ct) Al (p + a) Al (e, iti) 

12 



where B is algebraic on p. 

Finally, combining J\ with /, one gets 

■h = y/(B, Ri)(p~a)^(p + a)^(A,x) = 

(4.22) 

= (p - a) Xl/2 (p + a) Xl/2 (A,x) , 
- the first integral of ()3.6|) with A(p) algebraic on p. 

Obviously, J\ and J 2 are functionally independent. Indeed, the vectors A(p), B(p) are inde- 
pendent as being proportional to the lines £2,3 (p) of the matrix E _1 (p) whose determinant 
is not identically zero. Finally, in order to find the third linear first integral, we apply the 
Liouville theorem to the fundamental matrix solution of ([3.6)1 formed by the columns Xi(p), 
X2(p) and the arbitrary solution x(p): 

Z{p) = (X 1 {p),X 2 (p),x), 

(4.23) 

det(E) = (p- a) 2Xl (p + a) 2Xl det(Ri,R 2 ,x) = const • a(p) , 

where a(p) is a rational function of p in view of tr(A),tr(B) G Z. 
One derives from ()4.23|) the following first integral of (J3.6)) 

J 3 = (p - a)^ 2Al (p + a)- 2Al (C,x) , (4.24) 

with C(p) algebraic on p. 

In the case the vectors A,B,C are linearly independent, finding x from the linear system 
Ji = Ci, J 2 = c 2 , J3 = c 3 , ci >2 ,3 G C, we see that the general solution of (|3.6|) does not contain 
logarithmic branching and so we get a contradiction. Let us assume that A,B,C are linearly 
dependent. Then C = l\A + l 2 B where 1%, l 2 are certain algebraic functions of p and the 
following relation holds in view of (PT2T1) . (14^21 . (jQlj) 

(p-a)-i Al (p + a)^l Al ci/i + (p-ay Xl (p + a)- Xl c 2 l 2 = c 3 . (4.25) 

Since Ji, J 2 are functionally independent, the last expression shows that Ai G Q and the 
proof is finished with help of (j4.19j) . The case (A,Ri) 7^ 0, (B,R\) = is treated in the 
analogous way. □ 

With help of the Lemma El one obtains I(p, x + X 2 ) = J(p, x + (p — a) Xl (p + a) Xl R 2 ) which 
is a first integral of ([3.6)1 . As before we show (A,R 2 ) = (B,R 2 ) = 0. Since A, B are 
linearly independent, it follows from ()4.17j1 . 1)4.18)1 that X 2 = QX\ for a certain function 
6 = 6{p) ^ const. This situation can be ruled out by the substitution of Xi, X 2 into the 
system (J3.6)) that gives a contradiction. According to Lemma ^ the proof of Proposition El 
is finished. 

□ 
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5 Non-existence of additional analytic first integrals 

Our aim now is to prove the non-existence of new analytic first integrals of the rattleback 
problem (jl.7|) under one of the following hyperbolicity conditions 

Spectr(M±) = {s h2;3 G C : < |si| < 1, |s 2 | > 1, \s 3 \ > 1, s 2 ^ s 3 } , (5.1) 

or 

Spectr(M±) = {s 1)2 , 3 G C : > 1, < \s 2 \ < 1, < \s 3 \ < 1, s 2 ^ s 3 } , (5.2) 

which represent certain restrictions on the eigenvalues of the characteristic polynomial (|3.8|) . 
Recalling that £V A« = —1 and sis 2 s 3 = 1 we deduce from the relations e 2mXi = Si that at 
least one of the conditions ()5.1|) or ()5.2|) is satisfied when 

Ai,A 2 ,A 3 ^R and ImAi-ImA^O, Vi^j, (5.3) 

where Aj, i = 1,2,3 are three roots of the cubic algebraic equation 

x 3 + x 2 + d x x + fl = , (5.4) 

whose coefficients depend on m) and are defined by ()3.10|) . 

From now on, we will assume that the condition (j5.3J) is fulfilled and that the case (j5-H) 
holds. 

Let £ be the space of linear forms / = (L,x), L G C n dual to C n . To each M G GL(n, C) 
we associate the linear automorphism M : £ — > £ according to M J = (M T L, x). The next 
results shows how the hyperbolicity of G implies the reducibility of its polynomial invariants. 

Proposition 4. Lei the monodromy group G = (M + ,M_) have a polynomial homogeneous 
invariant P(xi, x, 2 , x 3 ). Then one of the following situations holds: 

a) P = x1x l 2 L(xi, x 2 , x 3 ), p, / G N, 

b) P = XiX^M(xi,x 2 , x 3 ), //, k G N, 

c) P = x\N(xi, x 2 , x 3 ), rj G N and M ± . Xi = SiXi, 
where L,M,N G C[xi, x 2 , x 3 ] . 

Proof. Since P is invariant under the action of M + and (|5.1|) holds, all different monomials 
entering in P contain a positive degree of x\. Therefore, P factorizes as follows 

P = x p 1 P 1 (x 1 ,x 2 ,x 3 ), p G N , (5.5) 

where P\ is a homogeneous polynomial not divisible by x\. 

Let £ be the set of all pairwise non-colliniear forms dividing P. In particular we have already 
X\ G S. We denote £q C CP 2 the set of directions of all elements from 8 equipped with 
the naturally defined action of G. The G-invariance of P implies clearly the invariance of 
directions from Sq under the action of G. 
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Exchanging the roles of M + and M_, one finds e £ £ - the eigenform of M_, M_ . e = Sie. 
Let e = eiXi + 62X2 + 63X3. Then, considering the orbit M™ . e, n > 1 and using (|5.1|) together 
with card^o < 00, one sees that if e\ 7^ then e2 = e3 = i.e. e = xi (case c). If e\ = 0, 
then either e 2 e 3 = and so e G {x 2 ,x 3 } (cases a-b) or e 2 e 3 7^ so that 3m G N such that 
S™ = s™. The last equation implies A 2 — A 3 G E which is impossible according to (|5.3|) . This 
concluded the proof. 

□ 

Proposition 5. The monodromy transformations M + and M_, taken in any basis, are 
permutationally conjugated i.e. 3C G GL(3,C) such that 

CM+C- 1 = M_, CM-C- 1 = M+ . (5.6) 

Proof. We take a basepoint e G 1Z on the positive imaginary axis Imp > 0. Let S(p), S(e) = 
Id be the normalised fundamental matrix solution of ()3.6|) and let G be the corresponding 
monodromy group. Since the equations (|3.6|) are invariant under the change of time p — > —p, 
S(p) = £(— p) is again a fundamental matrix solutions of (|3.6j) . Let r 1; T 2 be two loops 
starting from e and going around the singularities p = a and p = —a, respectively. We 
define the loops Tj = — 1\, i = 1,2 (symmetric to with respect to origin), starting from 
the point e = — e and having the same orientation as r^. Obviously, S(e) = Id and we can 
define the monodromy group G using E in the usual way. One sees that T\, T 2 define now 
the monodromy transformations around p = —a and p = a, respectively. The result follows 
then from the fact that G and G are always conjugated. □ 

Remark 3. It follows from the previous proposition that if C 2 7^ Id then it is a centraliser of 
G in GL(3, C). 

Proposition 6. Let us assume that the monodromy group G of the normal variational 
equations ()3.(i)l has a polynomial homogeneous invariant and that the conditions (jl.5j) . ()5.3j) 
hold. Then G is diagonalizable. 

Proof. One first assumes that a)-b) from Proposition 0] hold. 

Then we have three possible cases: £ = {xi,X2} (A), £ = {xi,X3} (B) or £ = {xi,X2,X3} 
(C). 

The group G acts on £ Q by permutations. Hence, exchanging if necessary X2 and X3, in A-B 
we can put 

M+ = diag(si, s 2 , s 3 ), M_ = 



at 
a 2 
a b cr 3 



, Oi G Spectr(M + ), a, b G C . (5.7) 



Indeed, £ always contains an eigenform of M_ corresponding to the stable eigenvalue si G 
Spectr(M + ) = Spectr(M_). In particular a\ = s\ or a 2 = s\. 

Using f)5.7j) one can calculate the unique (modulo multiplication by a constant) nonsingu- 
lar matrix T such that T" 1 M_ T = diag(si, s 2 , S3). Then, as follows from Proposition 
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T 1 M + T = M_. The last equation, together with ([5.3)1 and some elementary calculations 
implies a = b = that proves the result. 

If the case C holds, then the matrix M_ G SX(3,C) is either diagonal, so the proposition 
is proved, or is a permutation of the eigendirections of M + fixing one of them. In this last 
case, Spectr(M_) will contain necessarily a pair of eigenvalues with equal absolute values. 
This contradicts to ([5.3)1 . 

We consider now the case c) from Proposition 0] One puts P into the form 

P = x p 1 {x x Px{x l: x 2 , x 3 ) + P 2 (x 2 , x 3 )) , (5.8) 

where Pi, P 2 are homogeneous polynomials and P 2 ^ const since |si| < 1. 

Let M + = diag(s2, S3) and let M_ = (Tnij) 2 <i,j<3 denotes the restriction of the linear operator 
M_ to the x 2 , 23-plane. 

It is clear from ()5.8|) that P 2 is a polynomial semi-invariant for both M + and M + . If P 2 
contains two different monomials x^x™ and x^x^, n + m = p + q = degP — p then s r 2 = S3 
for r = n — p = q — m^O that contradicts to ([5.3)1 . 

If X2 (resp. X3) is the only monomial entering in P 2 , then x 2 (resp. x 3 ) is the eigenform of 
M_. 

Thus, exchanging if necessary x 2 and 23, it is sufficient to consider the case 



M + = diag(si,s 2 ,s 3 ), M_ 



sx 
a m 22 
b m 32 m 33 



(5.9) 



With help of ([5.6)1 . ()5.3[) and some elementary algebraic computations one shows that m 22 = 
s 2 , m 33 = s 3 . 

We introduce the matrices U = M+ X M_ and K = U — Id. One verifies that Spectr([7) = 
{1, 1, 1} and that if rang(f^) = 2 then M_ does not have any polynomial invariants. That 
can be done by transforming U to its Jordan form. If rang(i^) = then M_ = M + and 
the proposition is proved. The condition rang(J — Id) = 1 implies in turn: a = (i), 
b = m 32 = (ii) or m 32 = (iii). In the cases (i), (ii), M_ has the eigenform equal to x 2 
or x 3 that corresponds to the case considered before. If (iii) holds, it is sufficient to apply 
again the conjugacy conditions ()5.6[) to obtain a contradiction with (|5.3)1 . 

Finally, we consider the case then P 2 = cx^x 1 ^ t,l G N, c £ C*. Then M_ either preserves or 
permutes x 2 , 23-eigendirections of M + . In the first case the proceed as above. In the second 
one, it is sufficient to verify by a direct computation that Spectr(M_) will contain in this 
case a pair of eigenvalues ±s, s G C* so that the condition ()5.3[) is violated. The proof of 
Proposition |U1 is finished. 

□ 

The next proposition shows that in our case the Fuchsian system (13. 6 j) never has a diagonal 
monodromy group. 
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Proposition 7. Under the conditions (jl.5j) and (J5.3j) the monodromy group G of the normal 
variational equations ()3.6|) is not diagonalizable. 

Proof. Let us assume that G is diagonalizable: M + = diag(si, s 2 , s 3 ), M_ = diag(s ii; s i2 , SjJ, 
Sk = e 2mXk , k = 1, 2, 3 where (zi, z 2 , 23) is a certain permutation of (1, 2, 3). 

That implies existence of three independent solutions of (|3.6|) of the form 

X k (p) = (p - «) Afc (p + «) A ^iV fc (p), fc = 1, 2, 3 , (5.10) 
with vector functions Nk rational on p. 

Let Yi, i = 1, 2, 3 be three linearly independent eigenvectors of A with corresponding eigen- 
values A i (we remind that A, are pairwise different in view of (|5.3j) ) . One deduces from (|5.1U|) 
the following formulas containing the rational vector functions Rk 

X k = (p-a) x ^(p + a) x ^R k (p), k = 1,2,3 (5.11) 

with g k Uk G {1, 2, 3} and where now Rk = Y 9k + Y+(p — a) + ■ ■ ■ in the neighborhood of 
p = a and Rk = Y Uk + Y~{p + a) + ■ ■ ■ in the neighborhood of p = —a, (k = 1, 2, 3). Since 
the system (j3.6j) is invariant under the change p — > —p, Xk(—p) is also its solution. It yields, 
together with condition ()5.3|) . — Uk for k — 1, 2, 3. 

Under the same condition (j5.3J) . as a simple argument shows, Y gi , Y g2 , Y gA must be pairwise 
different. Since M + = M_ and M + M_M OQ = Id, the similar property holds for the point 
p = 00. 

Thus, one can represent Xk as below 

Xk = {p-a) Xk {p + a) Xk p n "Pk{p), k = 1,2,3, (5.12) 

where n k € Spectr(£>) = {0, —1, —1}; Pk are polynomial vector functions such that Pk(a) = 
Pk(—a) = Yk for k — 1, 2, 3 and -Pfc(O) are eigenvectors of 5. 

Let -Dfc be the order of Xk at infinity. One finds from ()5.12j) : Dk = — 2A& — nk — dk so that 
YjD k = - 2~2 n k - 2~2dk > -2 • (-1) - (0) ~Y.dk since Y n k < 0. Otherwise, 

we know that Y Dk = 4 (see Proposition |2J) and hence Yl dk < 2. Thus, at least one of 
the vectors Pk is constant and, as easy follows from substitution of ()5.12|) into ()3.6|) . is a 
common eigenvector of A and 5. This is impossible according to Proposition El and the proof 
is achieved. □ 

In view of Lemma ^ the above results can now be summarized by the main Theorem ^ . 



6 Conclusion 

The difficulty of the rattleback problem is due to its non-hamiltonian nature. That explains 
the non-trivial structure of the monodromy group studied in the previous sections. The tech- 
nical problem in applying our non-integrability Theorem are the hyperbolicity conditions 
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(J5.3j) . Of course, one can write these restrictions directly using Cardano formulas that will 
lead to quite complicated expressions. It may be interesting to consider a concrete example. 
Let Ji = 0.5, I 2 = 0.6, h = 0.8, 5 = 1.3, m = 1, g = 1, b x = 1, b 2 = 2, b 3 = 3. Then the 
eigenvalues of the characteristic polynomial ()3.8|) are 

\ x = -0.365 - 0.858i, A 2 = -0.435 + 0.963i, A 3 = -0.200 - 0.106i , 

and the conditions ()5.3|) are obviously satisfied. 

We believe that our non-integrability conditions can be strengthened. Thus, the further 
detailed analysis of the characteristic polynomial (|3.8|) is needed. 

We note that in the heavy rigid body case, the characteristic polynomial ()3.8)1 always has 
a real root as seen from (j3J). Indeed, it corresponds to existence of the fourth polynomial 
first integral (jl.llj) of the Euler-Poisson equations. In the rattleback case, the interesting 
remaining problem is the existence of new meromorphic first integrals. Our Theorem |21 
answers this question only then the variational equations ()3.6|) have logarithmic singularities. 
Our intention to avoid the study of the Zariski closure of the monodromy group G was 
twofold. Firstly, that makes the proofs self-contained and quite elementary. Secondly, we 
would like to underline the importance of the symmetry conditions (|5.6jh coming from the 
mechanical context of the problem and which simplify greatly the non-integrability analysis. 
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7 Appendix 

The matrix N is given by 



/ 



N 



i 2 ^ V 



p--a 
\ P 



v oy 



(7.1) 



The matrix M is given by 





( a iP 


a 2 p 






M = 




-aip 


C3 + C4 p2 _ Q 


(7.2) 




\ P 


p 2 —a 
P 


J 



The coefficients related by one quadratic equation, which is equivalent to ()3.14|) . 
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They are given by 

B d a x = 2 ^(pi - P2)(Si2 - hnpip 2 ) + (S12 - impip 2 )(Sn + E 22 + m(p\ - p\)) (7.3a) 

6 d a 2 = 2 "^"(Pi -P2)(S 22 + mpl) (7.3b) 

+ (E12 - iS 22 - impi(pi - p 2 )(Ei 2 + iE 22 + \mp 1 {p x - p 2 ) 

9 d a 3 = 2^(pi - p 2 )(Sn - mpl)- (7.3c) 

- (E n - iE 12 - mp 2 (pi - p 2 ))(E n + iS i2 + mp 2 {p 1 - p 2 )) 

26» d ci = 2m#(Ei 2 (pi - p 3 ) + iE 22 (p 2 - p 3 ) + impip 3 (pi - p 2 )) (7.3d) 

- m/3(Ei 2 (2pi - p 2 ) - iE 22 (pi - p 2 )) + im 2 fip x p 3 {p\ - p\) 

29 d c 3 = 2mg(T, 11 (p 1 - p 3 ) + iEi 2 (p 2 - p 3 ) - mp 2 p 3 (p x - p 2 )) (7.3e) 

- m/3(En(2pi - p 2 ) - iEi 2 (pi - p 2 )) - m 2 (3p 2 p 3 (p\ - p\) 

9 d c 2 = -m/3p 3 (Ei 2 pi + iE 22 p 2 ) (7.3f) 

6> d c 4 = m/3p 3 (E n pi + iEi 2 p 2 ) (7.3g) 



The matrix i? is given by 



The matrix A is given by 



-1 -2i(ci - c 2 ) 
5=10 -1 -2i(c 3 - c 4 ) 




-mi — ia 2 — ic 2 
-ia 3 iai — ic 4 
-i/p -1 



Rigid body case: 

ri = ~pi,r 2 = -ip 2 ,r 3 = p 3 = 0, m = 0,mg ^ 0, E 12 = S = . 
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